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I. INTRODUCTION 

Obtaining pure states through distillation processes is 
a key step to initialize and to control quantum systems 
in the field of quantum technology [T] [5]. For instance, 
state preparation is fundamental in several algorithms in 
quantum information and computation. So far various 
protocols for distillation have been proposed [21 [31 HI |S] . 

A procedure for distillation through Zcno-like measure- 
ments has been presented in [6] [J . The protocol there 
considered regards a generic bipartite system made of two 
interacting subsystems P and F. The unitary dynamics 
of the total system is interrupted at regular intervals by 
measurements performed on one of the two subsystems, 
say P. These measurements affect strongly the dynamics 
of the non-measured system F , which is then governed 
by some effective evolution operator. It is shown that if 
the spectrum of this operator satisfies certain conditions, 
then the non-measured system itself is driven toward a 
pure state irrespective of its initial conditions. This final 
state depends on the parameters of the total Hamilto- 
nian, the measurement one performs, i.e., the state on 
which the system is projected by the measurement, and 
the time interval between measurements. In this distilla- 
tion protocol measurements are performed repeatedly as 
in the case of quantum Zeno effect [8j [9l [TOl [TTl [T2] but 
here the interval between measurements is kept small but 
finite considering the limit of large number of measure- 
ments. 

This general procedure can be useful to initialize quan- 
tum systems and it has been also extended to the case 
where the system is made of many parts showing that 
it is possible to distillate entangled states [7l [T^. For 
example it has been shown that it is possible to estab- 
lish entanglement between two separate systems via re- 
peated measurements on a third entanglement mediator 
[14] . The effect of environmental noise on the protocol 
has been also investigated, showing that in the case of 
only dephasing effects one can still extract pure states 

However, these protocols have been analyzed only 
when the measured system has a discrete spectrum. Be- 



cause of the conditions required for distillation, it is not 
obvious that distillation can be obtained even when the 
measured system has a continuous spectrum. The main 
purpose of this paper is thus to apply the general proce- 
dure, as considered in [Bj, to the case where the measured 
system has a continuous spectrum. The system consid- 
ered is a particle (system P) , characterized by a continu- 
ous spectrum, interacting with a field mode (system F). 
Our protocol consists of measurements performed on the 
particle P to confirm it to be in a Gaussian state. We 
will show that the field mode is driven toward a pure 
state and we will examine how this distillation process 
depends on parameters like the mode frequency lo, the 
frequency of measurements r^^, the coupling between 
the particle and the field mode g, and the characteristics 
of the Gaussian state of the particle P to be confirmed 
repeatedly. 

The paper is organized as follows. In Sec. [IT] we de- 
scribe the features of the model we consider. In Sec. IIIII 



we describe the protocol and in Sec. IV we obtain and 
diagonalize the projected evolution operator. In Sec. [V] 
we discuss distillation in terms of physical parameters of 
the system. In Sec. [VI | w e summarize and discuss our 
results. In Appendices [A[ and [B] we collect some of the 
calculations, which are omitted in the text in order to 
keep the readability. 

II. MODEL 

We consider a particle of mass m interacting with a 
single field mode of frequency uj. The particle interacts 
linearly with the field mode and the Hamiltonian reads 
as 



2m 



+ hjj ( a) a + -z] + 9P{a^ + a) 



(1) 



where p is the particle momentum operator, a and a) 
are the annihilation and creation operators of the field 
mode satisfying the commutation rules, [a,a^] — 1, etc., 
and the real parameter g is the coupling constant. This 
Hamiltonian is analogous to the electromagnetic Hamil- 
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tonian restricted to the single-mode case and in dipole 
approximation |16} 117] . 

In the interaction picture the interaction Hamiltonian 
of Eq. ([T|) at time s takes the form 



H'{s) = gpiaU"^' + ae- 



(2) 



The evolution determined by this Hamiltonian can be 
treated exactly since its commutator at two different 
times, 

[H\s'),H^{s")] = -2i.9Vsincj(s' - s"), (3) 

commutes with the Hamiltonian itself. This allows one 
to obtain the exact evolution operator [THj as 



U'{t) ^ T^exp 



fi Jo 



dsH'{s) 



exp 



X exp 



2^ £ d. £ ds' 9{s - s')[H\s),H'{s')] 



^ Jo 



dsH^{s) 



(4) 



where is the time ordering operator and 0(s — s') is 
the Heaviside step function. Using Eqs. (|2|-([3| and 



ds / ds' 9{s — s') sin uj{s — s') 
Jo 



LOT — smwr 



, (5) 



the time evolution operator at time r can be put in the 
form 



IJ\t) = exp 



■ 2 -2 

ig p 



LOT — sm tjr 



X exp 



gp a 



it: 



- e^" 
Hlo 



Hlo 



(6) 



Indicating with Uq the time evolution operator associated 
to the free part of the Hamiltonian of Eq. ([T]) , the time 
evolution operator in the Schrodinger picture, U{t) — 
Uo{t)U'{t), is given by 



U{t) = exp 



i fP'T / 2mg^ lot — sinwr 
h 2m \ hco 



lot 



X exp 



X exp 



,„t" 1 
-iLOT I a' a + - 



9P a 



1 - e'" 

hco 



hio 



(7) 



III. PROTOCOL 

In this section we briefly review the distillation proce- 
dure based on repeated measurements as applied to 
our system. 

Assume that at time t = Q the particle is prepared in a 
pure state |<i>o) and the field in an arbitrary mixed state. 



denoted by p_f(0). The unitary dynamics of the total 
system, governed by the time-evolution operator U{t) in 
Eq. ([t]) , is interrupted by the measurements performed on 
the particle at intervals r. Each time, the measurement 
projects the particle in its initial state |<&o)- This action 
is represented by the projection operator O = \'^o){'^o\® 
Ip. The projection is partial on the total system because 
only the state of the particle is set back to its initial 
state while the field is not initialized, even though its 
dynamics is certainly affected by the measurements. The 
total system after N measurements is described by 

pTr{N) cx pU{T)f[\^^){'^^\®pF{mu\r)Of. (8) 

Following [6] , we introduce the projected evolution oper- 
ator between two consecutive measurements, 



Vr^{^o\U{r)\<^o 



(9) 



so that, after N measurements on the particle, the field 
is described by the density matrix 



N 



PF(o)i>r 

Pr{N) 



(10) 



Note that we retain only events in which the particle P 
is found in the state |<i>o) by every measurement. The 
normalization factor 



Pr{N) = TrpiV^^'pFm^''} 



(11) 



represents such a probability (up to N measurements) 
and gives the probability to obtain the state (10 1. 

In this paper, we consider the following protocol: the 
particle P is initially prepared in a Gaussian state |<i>o) 
and repeatedly projected on it at intervals r. Such a mea- 
surement would be realized by switching on a harmonic 
potential, whose ground state coincides with the Gaus- 
sian state |<i>o)j and seeing whether P is in the ground 
state. The Gaussian state |<i>o) is characterized by the 
variances of the coordinate Aro and of the momentum 
Apo, satisfying Arg Apo = ti/2. In the momentum space, 
it is given by 



I'&o)- / dp- 



1 



^M^poY 



exp 



P 



4(Apo)^ 



b), (12) 



where \p) are the eigenstates of the momentum operator 
p, the initial average momentum is po = 0, and the initial 
average position — Q. 

Once chosen the measurements to perform on the par- 
ticle, the next step is to ask if the protocol considered 
leads to purification of the field mode. In order to achieve 
purification for the non-measured system, the spectrum 
of the projected evolution operator Vr defined in Eq. (oj) 
must satisfy the following conditions [5]: its largest (m 
magnitude) eigenvalue must be unique, discrete, and non- 
degenerate. In the next section we shall compute the 
projected evolution operator Vr for the present setup and 
then diagonalize it in order to analyze its spectrum. 
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IV. PROJECTED EVOLUTION OPERATOR 



where 



Both the field evokition and the survival probability, 
given by Eqs. (10 1 and (11), depend essentially on the 
projected evolution operator Vr defined in Eq. (loj). This 
operator, using Eq. (12) as the state of the particle |4>o) 



to be measured repeatedly, results in 



Vr = exp 



'iLJT a' a 



1 



1 



dpe-^(f\ 
(13) 



where 



fip) = /3 



2(Apo)2 

\/2(l — cos Wt) ^ /, . N 

— -{a^e'^ + ae ^)p, (14a) 



+ / _ 2W^r-sinc.r \ ^^^^^ 

hm \ nuj LOT J 

Introducing three independent dimensionless parameters 



UJT, 



9 



and performing the integration in Eq. ( 13 1 we obtain the 
projected evolution operator as 



Vr = M exp 



X exp 



iuJT I a + 
2 ^ 



1 



ae-'^)2 



where 



M = 



l + zA2f [1 - 2g2 (1 _ 



G = 2AfVAp(l - cost) 



(16) 



(17a) 



(17b) 



A. Diagonalization 

As already stated, in order to check if the conditions for 
purification are satisfied in the present setup we need to 
analyze the spectrum of Vr- To diagonalize the projected 
evolution operator Vr it is useful to rewrite its expression 
in Eq. ( 16 1 in a single unified exponential. The details 
are given in Appendix wher e the projected evolution 
operator Vr is arranged in Eq. (AlOl in the form 



Vr = M CXp- 



ln(g - ^/q^ 



1) 



G 



5t2 



+ q [ a' a + 



(18) 



g = COST + sin T, g = cos r + — sin r. (19) 
It is diagonalized by the similarity transformation 

e^^e'^^'v^e-'^^'e-"^, (20) 



where 



C = g± - 1, ?7 = ±— , 

2q 



and then Vr is transformed to 
Vr M exp 



(21) 



(22) 



This shows that the eigenvalues 7„ and the right and left 
eigenstates, |u„) and (w„|, of 



Vr 



E 



ln\Un){Vn\ 



(15) are given by 



7„ = M exp 
\u„ 



T(n+l/2)ln(g- Vg2-1) 



e-^^'e-^^ln) 



(nle^V^', 



(23) 

(24a) 
(24b) 



where \n) is the eigenstate of the number operator d^a 
belonging to its eigenvalue n = 0, 1, . . . 

A comment is in order as to the choice of the sign 



(±) in Eqs. (21), (22), and (24). As shown in Appendix 
[B] in order to assure the normalizability of the eigen- 
states, for a given q, a correct sign must be chosen so 
that the inequahty \q ± y^q'^ — 1| = |g =F — 1|~"^ < 1 
holds. The real part of lii{q — q'^ — 1) is equal to 
ln|g — y/q"^ — l\ = — ln|g + y^^(f^^\ and it is easy to 
verify that one has to choose the upper(lower) sign, that 
is, +(-) in Eq. ^ and -(+) in Eqs. ^ and ^ if the 



real part is positive(negative). Therefore, ± may be sub- 
stituted with the sign of the real part of la{q— -y/g^ _ l). 
We observe that when f — {2£ + l)7r (^ = 0, 1, . . .) we 
have q = and then \q ± q^ — 1| = 1, independently 
of the choice of sign, so that for these values of r the 
normalizability of the eigenstates is lost and there is no 
distillation. 

V. DISTILLATION 

Analyzing the structure of the eigenvalues 7„ of Vr in 
Eq. ( 24 ) we see that the largest (in magnitude) eigenvalue 
7o is unique, discrete, and nondegenerate, so that in the 
large N limit, the operator Vj^ is dominated by a single 
term [6] 



V 



largo N 



7o l"o)(wo|- 



(25) 
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Then, using Eq. (25 1 in Eq. (10 1, in the large N hmit for 



a nonvanishing t, the state of the field asymptotically 
approaches the pure state, 



-rtAn large W Up Up ic\/cl 

Pp{N) * - — — - = 4 - 



(26) 



Notice that the pure state |uo) is explicitly written as 

\uq) ^e-'^^'^lO) = \/coshrS'(0|0) = VcoshrjO, (27) 

where S{^) = exp( — |a^'-^ — ^a?) is a squeezing operator 
with ^ = re*"^. The squeezing parameter r and the phase 
tp are given by 

r = tanh"^|C|, = arg C, (28) 



where C, is defined in Eq. (21 1 and q in Eq. (19 1. There 



fore, Eq. (26) shows that the field mode is distilled into a 
squeezed state. The final pure squeezed state |^) is inde- 
pendent of the choice of the initial state of the field, i.e., 
any initial (eventually mixed) state shall be driven to the 
unique pure state |^) by the repeated measurements per- 
formed on the particle to confirm it to be in the Gaussian 
state |<I>o). 

There are four, the first two independent of and the last 
two dependent on the initial state of the field mode, rel- 
evant quantities characterizing this distillation process: 

• Speed of distillation: the purification is achieved 
quickly if |7„/7o| = |7i/7or < 1 forn ^ 0. This 
means that we have a quick distillation if 



1 71/70 1 



exp 



Ink - 



1 



< 1. 



(29) 



The quickness of distillation is thus linked to 
In (7 



1 



ln|gf + \J - l\ 



• Degree of squeezing of the distilled state |^): this 
is given by r in Eq. (28). It regulates the av- 



erage number of quanta (a* a) = sinh r and its 
variance {{a)aY) — (a'a)^ = 2 sinh^ r (sinh^ r + 1). 
For a given r, (p determines the variance A? = 



{x^) of the quadrature operator xs 



^(ae + a^e**'), where -i? is a real phase. Ob- 
serve that e-272 < < e272. If si-a^i^ - <p) < 
(e^*" -I- 1)~^ the variance is less than 1/2 and 
in this case the quadrature x^ is squeezed. 

• Probability of success of the protocol: this is rep- 
resented by the survival probability Pt(N) intro- 
duced in Eq. (11). It is desirable to have a higher 



probability of success to keep higher yields when 
the distillation has been attained. While the state 



approaches the squeezed state |^) as shown in (26 1 
the probability behaves as 



PAN) |7oP^' 



uo\uo){va\pF{0)\va) 



(30) 



Therefore, it is preferable to have a larger I70I 
(closer to unity) for a slower decay of the proba- 
bility Pr{N). 



• Fidelity: it indicates how close the extracted state 
of the field mode is to the target pure state after N 
measurements. For an efficient distillation protocol 
the state after N measurements must be as close 
to the final target state as possible. The fidelity is 
defined by 



FAN) = muN)\o, 



(31) 



which is a positive number and approaches unity as 
the field state Pp{N) becomes close to the target 
pure state |^). 

In the following we shall analyze the dependence of 
these quantities on the parameters f, g, and Ap. Our 
aim is to find optimal values of the parameters satisfying 
two independent requirements: 

(i) fast distillation with high degree of squeezing [see 
(|29|], and 



(ii) high fidelity with a sufficiently high probability of 
success of the protocol [see (30)]. 



Even if the second condition is fundamental for an ef- 
ficient protocol for the distillation, we start discussing 
the first condition, concerning the speed of distillation 
and the degree of squeezing, because it involves quanti- 
ties which are independent of the initial state of the field 
mode. 

A. Distillation speed vs squeezing 

In order to have indications for the ranges of the pa- 
rameters where there is a quick distillation, we plot the 
ratio between the first two eigenvalues — In [71/70 1 = 

ln|g^ -^/g-^ — 1 1 , looking for regions where this quantity 

becomes large. Figure [l] shows — In [71/70 1 as a function 
of the dimensionless parameters in Eq. ( 15 ), f and g with 



fixed Ap — 1. The plot evidences that, for a fixed Ap, the 
ratio — In [71/70I has a strong dependence on the values 
of g and uj. In particular, for f < tt and g > l/\/2, the 
ratio can be greater than unity and in this region a fast 
distillation is available. 

Next , to investigate the degree of squeezing of the dis- 
tilled state, we plot the hyperbolic tangent of the squeez- 
ing parameter, tanhr — [C[, instead of the squeezing pa- 
rameter r itself, in order to avoid divergences in r when 
[CI 1. Figure [2] shows tanhr as a function of the di- 
mensionless parameters f and g with fixed Ap — 1. It 
is evident that for odd multiples of tt (values of r not 
allowed for distillation) and ^ 7^ we have tanh r ^ 1 
and thus r 00. With f just smaller than these values 
the distilled state is highly squeezed and in particular the 
region t < tt and g > l/-\/2 appears to be more appro- 
priate to obtain states with a high degree of squeezing. 
For g and Ap large enough, tanh r depends only on r. 

The above plots indicate that distillation speed, linked 
to — In [71/70 [, and the degree of squeezing, linked to 
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FIG. 1: The ratio between the first two eigenvalues 
— In 1 71/70 1 as a function of the dimensionless parameters f 
and g with fixed Ap = 1. 



tanh r 




FIG. 2: tanhr as a function of the dimensionless parameters 
f and g with fixed Ap = 1. 



tanh r, have strong and different dependencies on the 
parameters. Now we show that it is possible to find val- 
ues of the parameters where quick distillation and strong 
squeezing would be attainable simultaneously. In Fig. [3] 
we compare the behaviors of — In I71/70I (dotdashed line) 
and of tanhr (solid line), as functions of the dimension- 
less parameter f with g — I and Ap = 0.4. From the plot 
one sees that — In 1 71/70 1 > 0.5 (sufficiently quick distil- 
lation) and tanhr > 0.5 (sufficiently strong squeezing) 
can be fulfilled for tt/2 < lot < tt. 

B. Survival probability vs fidelity 

Both the survival probability and the fidelity depend 
on the initial state of the field. We consider the case 
where the field is initially in a coherent state \a). Al- 
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FIG. 3: — ln|7i/7o| (dotdashed line) and tanhr (solid line), 
as functions of the dimensionless parameter f with g — 1 and 
Ap = 0.4. 
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FIG. 4: Pt(N) (circles) vs Fr{N) (diamonds) as functions 
of the number of measurements A'' with g — 1, Ap = 0.4, 
and r = O.Qtt. The field is initially in a coherent state \a) of 
amplitude a — 1. 



though the survival probability and the fidelity can be 
analytically obtained when the field is initially in a co- 
herent state, their explicit expressions are rather involved 
so that in the following we will just present the results. 

Now, in order to check if the present protocol gives 
distillation with a good probability of success in this case, 
we compare the evolutions of the survival probability and 
of the fidelity for a given r. In particular we are looking 
for values of the parameters such that when the fidelity 
gets close to 1 the survival probability is still high enough. 
In Fig. |4] the evolutions of the survival probability and 
of the fidelity for a = 1 and f — O.Qtt are given. The 
values of the other parameters are the same as in Fig. [3] 
From the plot one sees that high values of fidelity can be 
obtained with the survival probability still far from 0. 

These two plots. Figs. |3] and |4] clearly indicate that our 
protocol allows one to generate efficiently, that is, with 
a high fidelity and a finite (nonvanishing) probability, 
pure squeezed states with a sufficiently high degree of 
squeezing. Indeed, for example, if we tune the parameters 
as ^ = 1, Ap — 0.4 and perform measurements with a 
period of f = O.Ott, we will obtain a well squeezed state 
with a squeezing parameter r ~ tanh~^ 0.6 ~ 0.54 with a 
fidelity of ^ 80% and a probability ~ 20% already after 
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a couple of measurements ~ 2, 3. 



APPENDIX A: UNIFICATION OF 
EXPONENTIAL FACTORS IN K 



VI. CONCLUSIONS 

In general for a bipartite system made of two inter- 
acting parts it is known that Zeno-like measurements 
on a part may lead, under certain conditions, the non- 
measured part towards a pure state independently of its 
initial configuration. The general procedure has been so 
far analyzed in the case where the measured part has a 
discrete spectrum, while it is not obvious whether the 
distillation can be obtained when the measured system 
has a continuous spectrum. Here we have analyzed this 
topic considering a specific bipartite system consisting of 
a particle, characterized by a continuous spectrum, inter- 
acting with a field mode characterized by a discrete spec- 
trum but with an infinite number of levels. The present 
distillation protocol consists of repeatedly projecting the 
particle to a Gaussian state. The projected evolution op- 
erator that regulates the field-mode dynamics between 
the consecutive measurements performed on the particle 
has been obtained. It has been shown that, with the 
measurement protocol chosen, the spectrum of this op- 
erator is discrete and satisfies the criteria that allow one 
to obtain a field-mode distillation. As a consequence of 
the protocol, the field is driven to a squeezed state in- 
dependently of its initial state. The dependencies of the 
distillation speed, that is connected to the ratio between 
the first two eigenvalues of the projected evolution oper- 
ator, and of the characteristics of the distilled state, i.e., 
the squeezing parameter, are investigated as functions 
of parameters such as the interval between two measure- 
ments, the particle-field mode coupling constant, and the 
width of the particle's Gaussian state. Varying the val- 
ues of the parameters different regimes are observed and 
we have shown that it is possible to choose values such 
that one has both quick distillation and strong squeez- 
ing and/or high values of fidelity and finite, well far from 
zero, values of probability of success of the measurement 
protocol. 
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The projected time-evolution operator Vr of Eq. (16 1 
is of the following form: 



where 



B = a)c 



(Al) 
(A2) 



1 

Note the Lee algebra among A, A\ and B, 

= [i,B] = 2i, [A\B] = -2A^. 

(A3) 

For these generators of the algebra, the following formula 
for the factorization of exponential is available: 



where 



< V = 



(/x/k) tanh k 
1 — (A/k) tanh k ' 

1 / sinh K 
— In cosh K — A 

2 V 

(y I k) tanh n 
1 — {\/ n) tanhAt' 



(A4) 



Inverse relations are given by 



(A5) 



/i = — xe 



sinh K 



sinh K 



ze 



K = cosh ^-(e^-|-e ^ — a;ze ^). 



A = - \J K? + [LV, 



Note the formula for the reciprocal function 
cosh^"'"a; = ± ln(a; + \J — 1). 



(A6) 



(A7) 



Now, let us come back to the projected time-evolution 
operator (All. By making use of the factorization for- 
mula (|A4|-(|A5[), one has 



% = Me-'"^-^e-TTcf<='""^^e-'"'^+'^)^e-TTG<="™"^ 

(A8) 

By exchanging the order of the first two exponentials, 
one obtains 



G 



Vr = Me- 



1 + G 



y = ln(l + G) + iujT. 



(A9) 



Then, we unify the exponentials via the formula (A4| 
with ( A6 1 to obtain 



K = Me'^(^'+-^)+^^, 
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G 



M = 



ln{q - ^q^ - 1) 



, G cos UJT + i sin ujt i-t. , 

A = ln((7- Vr - 1), 



(AlO) 



where q is defined in (19 1. The expression (18 1 is thus 
obtained. 

APPENDIX B: CHOICE OF SIGN IN 
DIAGONALIZATION 



Equations (211, (22 1, and (24 1 seem to show that the 



diagonahzation procedure presented in the text admits 
two possible signs ±. The immediate concern, in such 
a case, would be whether the unitarity of the time- 
evolution operator [/(r), which dictates that the abso- 
lute values of eigenvalues Y-r are strictly upper-bounded 
by unity, is preserved by the solution in Eq. (24 1, or 



which would be the right choice of the signs (Tjof the 
eigenvalues (and the eigenstates) if only one of them can 
be allowed. It will be shown here t hat th e normalizable 
eigenstates are those where \q — \/ q^ — V\ < 1, i.e., be- 
longing to the eigenvalues with their magnitudes always 
less than unity. 



doRdoi 



where the exponent is explicitly written as 



(B3) 



a ■ Aa — (ttR ai) 



1 + Cr, Ci 
Ci 1 - Cr 



Or 
ai 



(B4) 



and ttR — Re a, ai = Ima, Cr = ReC, Ci — ImC. The 
eigenvalues of the matrix A are easily found to be 1 ± 
both of which have to be positive, i.e., \(^\ < 1, in order 
for the above state is normalizable. That is, the nor- 
malizability of the eigenstates is assured if the condition 
Id < 1 is satisfied. This condition is explicitly written as 
that for q and G 



ICI 



1 <1, 



which just reduces to 



1 < 1. 



(B5) 



(B6) 



|7„|<exp {n+l/2)ln\q±y/q^ - l| 



< 1. 



(Bl) 



The normaliz ability of the eigenstate |uo) for n = 0, 

lle-^^'e-^-^lO)!! < oo, (B2) 

is sufficient to show the above statement. The left-hand 
side is calculated to be 



d a 



e 2 1." / 



la 



if G is replace with G^^. The normaliz ability condition 
of the eigenstates for the case of G^^ thus ensures the 
unitarity in the case of G. Stated differently, we have to 
choose an appropriate sign between -I- and — so that the 
absolute value of the argument of the logarithm satisfies 
the above inequality in order for the eigenstates to be 
normalizable. The unitarity is always satisfied, or we 
just have to make an appropriate choice of the phase of 
the square root ^/g^ — 1 = ^'(coswr -I- iG'sinwr)^ — 1, 
to which the sign (±) could be considered to be absorbed. 
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